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Understanding the oxidation of silicon has been critical to the success of all types of silicon materials, which 
are the cornerstones of modern silicon technologies. For the recent experimentally obtained 
two-dimensional silicene, oxidation should also be addressed to enable the development of silicene-based 
devices. Here we focus on silicene oxides (SOs) that result from the partial or full oxidation of silicene in the 
framework of density functional theory. It is found that the formation of SOs greatly depends on oxidation 
conditions, which concern the oxidizing agents of oxygen and hydroxyl. The honeycomb lattice of silicene 
may be preserved, distorted or destroyed after oxidation. The charge state of Si in partially oxidized silicene 
ranges from + 1 to +3, while that in fully oxidized silicene is +4. Metals, semimetals, semiconductors and 
insulators can all be found among the SOs, which show a wide spectrum of electronic structures. Our work 
indicates that the oxidation of silicene should be exquisitely controlled to obtain specific SOs with desired 
electronic properties. 

The progress of the experimental synthesis of silicene, the silicon counterpart of graphene, has been generat- 
ing great interest^"^. It is conjectured that silicene is even more advantageous than graphene^"^'*. For instance, 
the spin- orbit interaction of silicene is stronger than that of graphene, leading to detectable quantum spin- 
Hall effects in silicene^^'^^. The salient structural feature of buckling for silicene facilitates the tuning of the 
electronic properties of silicene by external fields^^"^^. In addition, silicene is clearly much more compatible with 
the existing well-established silicon technologies than graphene. Silicene-based devices may be readily incorpo- 
rated into current silicon integrated circuits. 

When silicene-based devices are fabricated, the oxidation of silicene is of the great question during the 
processing^^'^°. This is due to the fact that silicon routinely reacts with oxygen to form silicon oxide, which has 
been playing a critical role in the success of all types of current silicon materials such as bulk silicon^\ silicon 
nanowires^^'^^ and silicon nanocrystals^^'^^. It is rather reasonable that the oxidation of silicene should also be 
addressed to realize the full potential of silicene in all kinds of applications. De Padova et al.^^ and Molle et al.^° 
have recently shown that the oxidation of silicene only starts at a high dose (1000 L) of pure molecular oxygen. 
The relatively low reactivity of silicene to molecular oxygen may be due to the added sp^ hybridization of silicon 
atoms^^. When exposed to air, however, silicene is readily oxidized without identified reasons^°. This indicates 
that the mechanism for the oxidation of silicene may seriously depend on oxidation conditions. 

Up to now, very limited work has been carried out to investigate silicene oxides (SOs) that result from the 
oxidation of silicene, although it is imperative that the properties of SOs should be elucidated to guide the design 
of silicene-based device structures. By means of first-principles calculation, Wang et al.^^ have recently studied the 
mechanical and electronic properties of a specific SO with the stoichiometry of Si:0 = 1:1 (silicene monoxide). 
The silicene monoxide is obtained by analogizing its formation to that of a thermal-reduction-induced graphene 
monoxide^^. In fact, the graphene monoxide can only be produced by the thermal reduction of a multilayer 
graphene oxide. Diffusion in the multilayer structure limits the reduction process, leading to the kinetics-con- 
trolled formation of graphene monoxide^^. It is clear that the current non-existence of multilayer silicene oxides 
disables the formation of silicene monoxide in a fashion similar to graphene monoxide. In the present work, we 
begin with silicene, which is composed of a single layer of silicon atoms. Oxidizing agents such as atomic oxygen 
(O) and hydroxyl (OH) in a variety of bonding configurations are incorporated into the lattice of silicene to form 
SOs. The formation of all the SOs is evaluated in the point of view of thermodynamics. It turns out that the charge 
state of Si in partially oxidized silicene ranges from +1 to +3, while that of Si in fully oxidized silicene is +4. 
When silicene is oxidized by O, single-side over-bridging O (O^^) is the most likely incorporated O configuration 
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during the oxidation. If OH is the oxidizing agent, boat-Uke OH 
(OHhi) and umbrella-Uke OH (OH^/) may be the most hkely incor- 
porated in the resulting SOs. When O and OH both exist, fully 
oxidized silicene may be readily produced. The structural character- 
istics of all the SOs have been identified in this work. We find that the 
electronic properties of SOs significantly depend on the bonding of O 
and OH. Metallic, semimetallic, semiconducting and insulating SOs 
can all be obtained. 

Results 

Partially oxidized silicene. Figure 1 shows the optimized structures 
of silicene and SOs that result from the partial oxidation of silicene 
with both top and side views. The characteristics of all the optimized 
structures are tabulated in Table 1 (notations of all the O and OH 
configurations are introduced in Table SI in the Supporting 
Information). In the single-atom-thick honeycomb lattice of 
silicene, the bond length of Si-Si (2.28 A), the bond angle of Si-Si- 
Si (116°) and the buckling distance (0.45 A) are all consistent with 
those reported in previous works^^'^^"^° (Figure 1 (a)). When each Si- 
Si bond is transformed to a Si-O-Si bond by an O atom during the 
oxidation of silicene, the buckling of Si atoms actually disappears in 
the resulting SO (Figure 1 (b)). All the O atoms are also in the same 
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Figure 1 | Optimized structures of (a) silicene and silicene oxides with (b) 
Oqb. (c) Odb, (d) of,, (e) O^,,, (f) OH,;, (g) OH^, (h) OH^and (i) OH^. 
Both top and side views of the structures are shown. The unit cell of each 
structure is indicated by a rhombus or rectangle. Si, O and H atoms are 
denoted by green, red and grey balls, respectively. 



plane, which is only 0.02 A away from that of Si atoms. Therefore, 
quasi-in-plane bridging O (O^^) can be used to denote the 
incorporated O. It is likely that each Si-Si bond at only one pair of 
opposite sides in a hexagonal ring is oxidized. But the oxidation 
results in two Si-O-Si bonds, which are symmetric with respect to 
the oxidation-flatted plane of Si atoms (Figure 1 (c)). This leads to the 
configuration of double bridging O (O^^) for the incorporated O. An 
O atom may also overbridge two neighboring Si atoms via their 
unpaired electrons. The overbridging O atoms can be located at 
either both sides (Of^) or single side (OJ^) of the original silicene 
(Figure 1 (d) and Figure 1 (e)). When the oxidizing agent is OH, each 
Si atom may be simply bonded to an OH group via the unpaired 
electron of Si. For a single hexagonal ring in the honeycomb lattice of 
silicene, the number of OH groups located at one side of the 
hexagonal ring may be zero, one, two or three. We denote the OH 
groups in the resulting SOs as top-like OH (OH^/) (Figure 1 (f)), 
umbrella-like OH (OH^/) (Figure 1 (g)), boat-like OH (OH^/) 
(Figure 1 (h)) and armchair-like OH (OH^/) (Figure 1 (i)), 
respectively. Please note that OH^i and OH^/ alternately passivate 
Si atoms with a step of one and two Si atoms at both sides of the 
original silicene, respectively. Among all the oxidizing agents only 
OHai and OH^/ do not change the original arrangement of Si atoms in 
silicene. The buckling distance, Si-Si bond length and Si-Si-Si bond 
angel of the SO with OH^/ or OH^/ are all the same as those of silicene 
(Table 1). In the lattices of the SOs with OH^/ and OH„/ all the Si-Si-Si 
bond angles decrease to be between 103 and 113°. The Si-Si bonds in 
the SO with OH^/ are either stretched or compressed, while those in 
the SO with OH„/ are all stretched. For the SO with O^p the bond 
length of Si-Si and bond angle of Si-Si-Si decrease to 2.19 A and 106°, 
respectively. All the Si-Si bonds are stretched and some of Si-Si-Si 
bond angles increase up to 122° in the SO with OJ^ or O^^. All the 
above-mentioned structural changes lead to the distortion of the 
honeycomb lattice of silicene. The buckling distances of Si atoms 
in the distorted honeycomb lattices of the SOs with OH^/, OH^/, 
Odh OJ^ and Of^ are 1.14, 0.03, 0, 0.02 and 0.03 A, respectively. It 
is clear that no Si-Si bonds exist in the SO with O^^ (Figure 1 (b)). The 
length of —83% Si-O bonds is 1.91 A, while that of the remaining Si- 
O bonds is 1.42 A. The Si-O-Si bond angles fall in the range from 104 
to 112°. All the structural characteristics of the SO with O^^ indicate 
that the honeycomb lattice of silicene is destroyed by the incorpo- 
ration of Oqh. 

We have worked out the formula of each SO by examining its unit 
cell (Table 1). The formulas of the SOs with O^^, Odb, O^^, O^^, OH,/, 
OH^/, OHyi and OH^/ are SisOj, Si404, SigOj, SigOj, Si2(OH)2, 
Si6(OH)6, Si4(OH)4 and Si2(OH)2, respectively. Therefore, the charge 
states of Si are +3 and +2 for the SOs with O^/, and O^by respectively. 
The charge states of Si are all +1 for the SOs with O^^, O'^^, OH,/, 
OHuh OHbi and OH^/. It is clear that the highest charge state (+4) of 
Si has not been obtained in these SOs. Therefore, they only result 
from the partial oxidation of silicene. 

The formation energy (Ey) per atom of oxygen in a SO can be 
calculated by using 



No 



Ef=[Et 



where E, is the total energy of the SO. E^ is the total energy of silicene. 
No and Nqh are the numbers of O and OH, respectively. fiQ^ and 
are the chemical potentials of O2 and OH, respectively. Figure 2 
shows the change of Ey with respect to A^q^ for each SO by fixing 
at the total energy of oxygen gas ( — 7.5 eV). A/Hq^ is the difference 
between fiQ^ in the oxidation system of silicene and that in water 
(—10.8 eV). For the SOs induced by the partial oxidation with O (e. 
g., in the atmosphere of oxygen gas), Ey increases in the order of Of^, 
O^^, Odb and O^b- For those induced by the partial oxidation with OH 
(e. g., in water), Ey increases in the order of OH/,//OH„/, OH^i and 
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Table 1 | Formulas and structural properties such as the buckling 
distance (A) of Si, the bond length of Si-Si and the bond angle of Si- 
Si-Si forsiliceneand partially oxidized silicene(SOswith Oqb, Ojb, 
Oi' OH,/, OH,/, Ohu and OH,/). Each bandgap (Eg) is 
obtained with the PBE correlation exchange functional. The more 
accurate values of Eg obtained with the B3LYP correlation 
exchange functional are shown in the brackets 





Formula 


A (A) 


<Si-Si> (A) 


zi(Si-Si-Si) (°) 


^g(eV) 


Silicene 


Si2 


0.45 


2.28 


116 


0(0) 




SisOs 


0 






5.2 (6.6) 


Odb 


Si404 


0 


2.19 


106 




Of, 


Si603 


0.03 


2.34 


1 10-122 


0.5 (0.7) 


o:, 


Si603 


0.02 


2.34 


1 10-122 


0.3 (0.6) 


OH^ 


Si2(OH)2 


0.45 


2.28 


116 




OH„, 


Si6(OH)6 


0.03 


2.34 


104-1 13 


0.2 (0.4) 


OHw 


Si4(OH)4 


1.14 


2.06-2.34 


103-1 12 


0(0) 


OH„, 


Si2(OH)2 


0.45 


2.28 


116 


0.4 (0.4) 



OHti. Therefore, we should expect that the oxidation of siUcene in the 
atmosphere of oxygen gas the most hkely leads to the formation of 
the SO with O^^. When silicene is only oxidized by OH, the SOs with 
OHyi and OH^/ are the most likely produced. 

We now move to the electronic properties of all the partially oxi- 
dized silicene. Figure 3 shows the band structure and density of states 
(DOS) of each partially oxidized siHcene. For comparison, the band 
structure and DOS of silicene have also been included in Figure 3. 
Consistent with previous results, silicene is semimetallic with a zero 
bandgap^ ^'^^'^^ (Figure 3 (a)). The top of the valence band and the 
bottom of the conduction band linearly cross at the Fermi level, 
where the DOS is negligible. O^^ opens the bandgap of silicene by 
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Figure 2 | Formation energy of a silicene oxide with O^f,, Ojf,, O^^, O^^, 

OH,;, OH„;, OH^;, OH,;, O^;, + Oj;,, O^;, + O^,^, O^;, + Of^, O^^ + OH,;, 

Oqb + OHu;, Oqiy + OH;,; or O^;, + OH^ with respect to the difference 
between the chemical potential of OH in the oxidation system of silicene 
and that in water. 

5.2 eV (Figure 3 (b)). Such a large bandgap means that the SO with 
O^^ is actually an insulator. We notice that there are two deep energy 
levels in the bandgap of the SO with O^h- By examining the partial 
DOS of the SO with O^^ (Figure SI in the Supporting Information), 
we find that the two deep energy levels mainly originate from the p 
orbital of O^^. Both the p-orbital of Si and the p-orbital of O signifi- 
cantly contribute to the conduction band of the SO with O^^. The 




Figure 3 | Band structures and density of states (DOS) of (a) silicene and a silicene oxide with (b) O^^, (c) Oj^, (d) O^^, (e) O^^,, (f) OH,;, (g) OH^;, (h) 
OH;,; or (i) OHa;. F, K, M are special points in the first Brillouin zone of a cubic crystal. F, X, S are special points in the first Brillouin zone of a 
hexagonal crystal. Energy is shifted so that the Fermi level is 0 eV (the horizontal dashed line). Conduction bands, valance bands, deep energy levels and 
DOS are indicated by red, blue, green and pink lines, respectively. 
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valence band of the SO with O^^ mainly originates from the p-orbital 
of O. It is known that a bandgap is usually underestimated in the 
density functional theory (DFT) based on the Perdew-Burke- 
Ernzerhof (PBE) correlation exchange functional. Thus, we have also 
performed DFT calculations based on the Becke, three-parameter, 
Lee-Yang-Parr (B3LYP) exchange- correlation functional to get the 
more accurate bandgap. For the SO with O^^, the B3LYP bandgap is 
6.6 eV (Table 1). 

The Fermi level is located in the conduction band for the SO with 
(Figure 3 (c)), indicating that the SO with O^^ is a metal. Figure 3 
(d) shows that the SO with Of^ is an indirect-bandgap semi- 
conductor. The valence -band maximum appears at T point, while 
the conduction-band minimum appears at M point. The bandgap is 
as large as 0.5 eV (the B3LYP bandgap is 0.7 eV). The SO with O^y is 
also an indirect-bandgap semiconductor (Figure 3 (e)). The valence- 
band maximum and conduction-band minimum are located at T 
point and K point, respectively. The bandgap is 0.3 eV (the B3LYP 
bandgap is 0.6 eV). There are no energy levels in the bandgap of the 
SO with Of^ or O^^. The origins of the conduction band and valence 
band for the SO with O^^, Of^ or O^^ are basically similar to those for 
the SO with O^^ (Figure SI in the Supporting Information). 

For the SO with OH^/ the Fermi level is located at the bottom of the 
conduction band (Figure 3 (f)). The DOS near the Fermi level is 
relatively large. Therefore, the SO with OH^/ is actually a metal. For 
the SO with OH„/, the valence-band maximum and conduction-band 
minimum are located at K' point (between K point and T point) and 
M point, respectively (Figure 3 (g)). The bandgap is 0.2 eV (the 
B3LYP bandgap is 0.3 eV). As shown in Figure 3 (h), the Fermi level 
crosses both the valence-band maximum and conduction-band 
minimum for the SO with OH^/. The DOS at the Fermi level is 
negligible. Clearly, the SO with OH^/ is also a semimetal. The val- 
ence-band maximum and the conduction-band minimum of the SO 
with OHai are located at M' point (between T point and M point) and 
r' point (between T point and M point), respectively (Figure 3 (i)). 
The indirect bandgap of the SO with OH^/ is 0.4 eV (the B3LYP 
bandgap is also 0.4 eV). The valence band of the SO with one type 
of all the OH configurations also mainly originates from the p-orbital 
of O. The s-orbital of H now significantly contributes to the conduc- 
tion band together with the p-orbital of O and the p-orbital of Si 
(Figure SI in the Supporting Information). 

For all the semiconducting partially oxidized silicene (SOs with 
O^^, O^^, OHui and OH^/), we have calculated the effective masses 
(m*) of electrons and holes by using the method adopted by Ni et al.^^ 
and Li et al.^^ (Table S2 in the Supporting Information). It is found 
that m* is between 0.36 and 4.78 mo (mo is the electron rest mass). 
According to /i = ei/m* (^ is the carrier mobility, e is the charge of an 
electron, t is the scattering time that is assumed to approximate to 
that in graphene), we obtain that the carrier mobilities of the semi- 
conducting partially oxidized silicene are in the order of magnitude 
of 10-10^ cm^ V"^ s"^ (Table S2 in the Supporting Information). 

Fully oxidized silicene. We find that O^h should be incorporated 
together with other types of O or OH configurations to achieve the 
charge state of + 4 for Si in fully oxidized silicene. Figure 4 shows the 
SOs based on the combination of O^^ and one of the configurations 
of Odh Of^, Oly, OHti, OHuh OHu and OH«/. Please note that the 
number of O^^ atoms is reduced by one third for the SO with Oqy + 
Odb because the existence of O^b only allows four Oqb atoms to be 
incorporated in a hexagonal ring. The buckling distances of Si atoms 
in the SO with O^b + O^b and O^b + Of^ are 0.04 and 0.06 A, 
respectively. Si atoms in all the other SOs are located in the same 
plane. This is the same as what happens to the SO only with O^b- By 
examining the unit cells (Figure 4), we obtain that the formulas of the 
SOs with Oqb + Odb. Oqb + Of^, Oqb + O^^^, Oqb + OH,,, O^b + OH,,, 
Oqb + OH,,, and O^b + OH^, are Si204, SigO^, SigO^, Si406(OH)4, 



(a) 0,, + 0^^ (b) + (c) O.t'^O'ot 




(d)0^^ + OH,/ (e) O^^ + OHu, (f) O^^ + OH^/ (g) O^^ + OH^/ 




Figure 4 | Optimized structures of silicene oxides with (a) Ogf,+ Ojf„ (b) 
Oqb + Of^, (c) + O^,^,, (d) + OH,;, (e) O^;, + OH„;, (f) O^;, + OH^ 
and (g) Ogb+OHa/. Both top and side views are shown. The unit cell of 
each structure is indicated by a rhombus or rectangle. Si, O and H atoms 
are denoted by green, red and grey balls, respectively. 

Si609(OH)6, Si203(OH)2 and Si203(OH)2, respectively (Table S3 in 
the Supporting Information). 

The formation energy of each fully oxidized silicene is also shown 
in Figure 2. When silicene is only oxidized in the atmosphere of 
oxygen gas, the most likely formed fully oxidized silicene is the SO 
with Oqb + O'^ij. If silicene is oxidized in an environment in which 
both O and OH exist, a high chemical potential of OH may most 
likely lead to the formation of the SO with Oqb + OH„. 

Figure 5 shows the band structure and DOS of each fully oxidized 
silicene. For the SO with Oqb + OH,,, the energy levels introduced by 
OH,, can basically fill up the bandgap that is opened by Oqb. 
Therefore, the full oxidation of silicene with Oqb + OH,, gives rise 
to a metalUc product. The PBE (B3LYP) bandgaps of the SOs with 
Oqb + Odb. Oqb + OH^, and Oqb + OHai are 5.1 (6.4), 5.0 (5.4) and 
5.1 (7.2) eV, respectively. This indicates that the SOs with Oqb + O^h 
Oqb + OH,^, and Oqb + OH^, are insulators. Nevertheless, we should 
notice that there are deep energy levels in the rather wide bandgaps of 
these three insulators (Figure 5 (a), (f) and (g)). The SOs with Oqb + 
O^^, Oqb + O^^ and Oqb + OH„ are all indirect-bandgap semi- 
conductors. Their PBE (B3LYP) bandgaps are 2.8 (3.9), 3.2 (4.3) 
and 2.1 (3.6) eV, respectively. Again, we see deep energy levels in 
the bandgaps of these three semiconductors. It should be pointed out 
that both the valance bands and conduction bands of all the fully 
oxidized silicene mainly originate from the p-orbital of O (Figure SI 
in the Supporting Information). This signifies that O plays the dom- 
inant role for the electronic properties of all the fully oxidized sili- 
cene. The effective masses of electrons and holes for the 
semiconducting fully oxidized silicene are between 0.73 and 87.81 
mo, leading to carrier mobilities in the order of magnitude of a few to 
10^ cm^ s~^ (Table S2 in the Supporting Information). 

Discussion 

De Padova et al.^^ and Molle et al.^° have recently carried out experi- 
ments on the oxidation of silicene in the atmosphere of oxygen gas. 
The low charge state (+1 and +2) of Si observed after the high 
oxygen exposure (1000 L) of silicene may be attributed to the forma- 
tion of the SOs with O^^, Of^ and O^h which have lower Ef. The Eyof 
the SO with Oqb is even higher than those of the SOs with Oqb + O^^ 
and Oqb + O^^, and comparable to that of the SO with Oqb + O^b- 
This means that +3 cannot be the dominant charge state for Si in 
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Figure 5 | Band structure and density of states (DOS) of a silicene oxide with (a) Oqb+Odb, (b) Oqb + O^y, (c) Oqb + O^^, (d) Oqb + OHfi, (e) 0^^, + 
OHu7, (f) Oqf, + OHbi or (g) Ogf, + OHg;. F, K, K', M are special points in the first Brillouin zone of a cubic crystal. F, X, S are special points in the 
first Brillouin zone of a hexagonal crystal. Energy is shifted so that the Fermi level is 0 eV (horizontal dashed line). Conduction bands, valance bands, deep 
energy levels and DOS are indicated by red, blue, green and pink lines, respectively. 



only O-oxidized silicene. The charge state of Si may directly change 
from + 1/+2 to +4 with the advancement of silicene oxidation in the 
atmosphere of oxygen gas^^. Please note that silicene cannot be fully 
oxidized if it is only oxidized by OH (Figure 2 and Table 1). The 
charge state of Si in the final OH-induced SO is always + 1 unless H is 
detached from OH to give rise to restructuring. 

When silicene is oxidized in air, O from oxygen molecules and OH 
from water may both account for the oxidation. It appears that in 
MoUe et al.'s work^° the air was an OH-rich environment, which 
lowered the Eyof each OH-incorporated SO. The oxidation of sihcene 
was easily initiated to produce the SOs with all types of OH config- 
urations, in which the charge state of Si was all + 1 . O^^ then joined in 
the oxidation, leading to fully oxidized silicene by means of the 
combination of O^^ and one type of OH configuration. Therefore, 
the charge state of Si changed from + 1 to +4 during the oxidation of 
silicene in air^°. 

We have shown that a variety of SOs may be produced by the 
oxidation of sihcene. The spectrum of electronic properties assoc- 
iated with these SOs is even wider than that obtained by the oxidation 
of graphene^^. On one hand, we have been provided with an oppor- 
tunity of making vastly different materials such as metals, semime- 
tals, semiconductors and insulators from the single material of 
silicene only by means of oxidation. On the other hand, excellent 
control on the oxidation conditions of silicene needs to be exerted to 
obtain specific SOs with desired electronic properties. For example, 
one can control the oxidation of silicene in the atmosphere of oxygen 
gas to produce partially oxidized silicene with O^^ or Of^. This effec- 
tively opens the bandgap of silicene without introducing any deep 
energy levels in the bandgap (Figure 3). The resulting semiconductor 
(the SO with O'^^ or Of^) with a bandgap of 0.6 or 0.7 eV can very well 
remedy the zero-bandgap disadvantage of silicene. We should note 
that this bandgap opening is much more significant than those 



induced by vertical electric fields^^'^^. By exquisitely tuning the oxida- 
tion of silicene with both O and OH, one may also produce insulators 
such as the SOs with O^^ + OH^/ and O^^ + OH^/. These insulating 
SOs can be readily used in silicene-based devices such as field-effect 
transistors as gate oxides^^'^l 

Methods 

Density functional theory (DFT) as implemented in the Vienna ab initio simulation 
package (VASP) with the projector-augmented wave (PAW) method^^'^^ has been 
employed in this work. The Perdew-Burke-Ernzerhof (PBE) correlation exchange 
functional at the generalized gradient approximation (GGA) level is adopted. The 
structures of silicene and SOs are initially optimized at 0 K. During the optimization 
the total energy and atomic forces for each structure are minimized. The sampling of 
reciprocal space is performed with the Monkhorst-Pack scheme^^. The meshes of (5 X 
5X1) and (9 X 9 X 1) special k points are used for the structural relaxation and 
density- of- states (DOS) calculation, respectively. Depending on the symmetry of the 
unit cell for a specific structure, an appropriate path in the first Brillouin zone is 
adopted to calculate the band structure^^. All the calculations proceed until the 
changes in energy and the force on each atom are less than 1 X 10"^ eV per cell and 1 
X 10"^ eV/A, respectively. The Becke, three-parameter, Lee-Yang-Parr (B3LYP) 
exchange-correlation functional has also been used to calculate the bandgaps of SOs 
after their structures are optimized. 
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